Single qubit gates with Jump and Return sequences 
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We discuss the implementation of frequency selective rotations using sequences of hard pulses 
and delays. These rotations are suitable for implementing single qubit gates in Nuclear Magnetic 
Resonance (NMR) quantum computers, but can also be used in other related implementations 
of quantum computing. We also derive methods for implementing hard pulses in the presence 
of moderate off-resonance effects, and describe a simple procedure for implementing a hard 180° 
rotation in a two spin system. Finally we show how these two approaches can be combined to 
produce more accurate frequency selective rotations. 

PACS numbers: 03.67.Lx,82.56.-b 



I. INTRODUCTION 

A key element in any implementation of quantum com- 
puting is the ability to perform single qubit gates. 
Clearly these require the ability to rotate one qubit (spin) 
while leaving the rest of the system unaffected. In many 
proposed implementations this is achieved using the fact 
that qubits are spatially localised, but in liquid state Nu- 
clear Magnetic Resonance (NMR) systems 2j this is not 
possible, and it is instead necessary to distinguish the 
qubits by their different resonance frequencies. Thus fre- 
quency selective excitation plays a key role in NMR quan- 
tum information processing 0, 0, IE IE Hi ■ 

The simplest form of excitation in an NMR experiment 
is a short high power burst of radiofrequency radiation, 
known as a "hard" rf pulse. Such a pulse will excite spins 
over a considerable range of frequencies: a simple (if not 
entirely correct) argument is that the frequency do- 
main excitation profile is the Fourier transform of the 
time domain pulse shape. Frequency selective excitation 
could in principle be achieved by using a long low power 
pulse (a "soft" pulse), but the selectivity is poor as the 
Fourier transform of a rectangular pulse is a sine func- 
tion, which has considerable power in the wings. A better 
method is to use shaped pulses, instead of constant power 
rectangular pulses, thus producing a more desirable pro- 
file in the frequency domain . 

This approach has been extensively developed, both for 
conventional NMR experiments (9l and for NMR quan- 
tum information processing 0, [n| (note that the prop- 
erties of many modern shaped pulses cannot be rational- 
ized using Fourier arguments and it is frequently neces- 
sary to use numerical methods to calculate them). In 
some experiments, however, it is overcomplex. Consider 
the case where a hard excitation pulse will only affect 
two qubits, for example in an NMR system which only 
has two spins of a particular nuclear species. In this case 



a much simpler scheme based on hard pulses and delays 
can be applied [HllEIIl- 



II. JUMP AND RETURN 

Schemes for frequency selective rotation using hard 
pulses are ultimately based on the "Jump and Return" 
solvent suppression pulse sequence ,15^ JJJ , which is de- 
signed to avoid exciting signals directly resonant with 
the applied radiation while (at least partially) exciting 
all other spins. To see how the Jump and Return se- 
quence can be modified for selective rotation, it is helpful 
to show how this sequence works for solvent suppression, 
which is effectively selective non-rotation. 

The Jump and Return sequence comprises a period of 
free precession sandwiched by two 90° hard pulses on 
resonance with the solvent frequency. 
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This and subsequent pulse sequences are written in NMR 
shorthand notation, where 90 is the nutation angle in de- 
grees and subscripts indicate the nutation axis, which for 
ideal pulses lies in the xy-plane; note that time runs from 
left to right, so that the 90_^ pulse is applied first. We 
assume that the system begins in its thermal equilibrium 
state, with all the spins along the z-axis of the Bloch 
sphere 9]. The first hard pulse excites all the spins, and 
during the period of free precession any spins away from 
the solvent frequency will precess. The solvent spins do 
not precess in the rotating frame 0, |E| , and are returned 
to the initial state by the second pulse, but non-solvent 
spins that have strayed during the period of free pre- 
cession are not fully returned, and hence result in some 
excitation (see figure ^1. It is easily shown that the se- 
quence in Eq. ^is equivalent to a O^y hard pulse, where 
is the precession angle during the period t. 

The amount of excitation that a spin receives thus de- 
pends on the offset of the spin frequency from the pulse 
frequency, as for a given duration between pulses, the ex- 
tent of the offset corresponds to an angle of rotation in 
the xy-pXauB. A spin will be maximally excited if it has 
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rotation. In particular 



90_ 



FIG. 1: Jump and Return solvent suppression sequence: the 
solvent signal (at zero frequency in the rotating frame) expe- 
riences no overall evolution, while spins at other frequencies 
are partially excited. 
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where we use non-selective y-rotations (hard pulses), to 
convert a contra-axial z-rotation (an appropriate period 
of free precession) into a contra-axial x-rotation. Rota- 
tions with other phase angles can, of course, be achieved 
by changing the phases of the two outer pulses. For a 
selective, rather than contra-axial, xy-rotation we can 
apply another non-selective a;?/-rotation along the same 
axis as the contra-axial xy-rotation |12| | , so that 



precessed 90° from its initial position in the xy-j>\a,ne 
(and therefore is completely unaffected by the second 
pulse), as shown in figure ^ There are two offsets where 
spins will precess 90° in the xy-plane {i.e., be maximally 
excited) and these are equally distant above and below 
the pulse frequency. Spins at these positive and negative 
offsets will end up along the same axis in the xy-plane but 
pointing in opposite directions, corresponding to overall 
positive and negative rotations respectively. 

Now we can return to our original situation with only 
two spins and consider how the Jump and Return se- 
quence can be useful for selective rotation. The reference 
frequency is set halfway between the resonance frequen- 
cies of the two spins, so that in the rotating frame the 
spins will precess at equal but opposite rates, ±(5tJ. The 
Hamiltonian can then be written using product operator 
notation ^3 



(2) 



where J is the spin-spin coupling and we have assumed 
a weakly coupled system, where \5u)\ » \J\. The effect 
of weak coupling is to replace each transition frequency 
by a cluster of transitions, known as a multiplet, but 
if the width of the multiplet is narrow compared with 
the separation between multiplets (which will be true for 
weak coupling) then the effects of spin-spin coupling can 
be neglected. We also assume for the moment that the 
hard pulses are instantaneous, so that there is no pre- 
cession during a pulse (we will return to this assumption 
later). If we set the delay period between the pulses to 
r = 7r/25a; the two spins will rotate by ±90° during the 
precession period, corresponding to the offset resulting 
in maximum excitation in the solvent suppression case. 

So far we have largely considered excitation sequences, 
that is sequences designed to take a spin from the z- 
axis into the xy-plane, but these sequences can be used 
to perform the general rotations required for quantum 
computing, as the Jump and Return sequence can be 
viewed as a simple rearrangement of the composite z- 
pulses (sequences designed to perform z-rotations using 
only xy-rotations) devised by Freeman et al. In 
general, one can perform any rotation around the Bloch 
sphere using rotations that are orthogonal to the desired 
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where 1 indicates the identity (do-nothing) operation. 
The extra non-selective rotation will return the spin we 
do not wish to rotate to its initial position, and further 
rotate the spin we want to select. Note that we halve the 
initial contra-axial rotation angle to compensate for the 
extra rotation on the selected spin. 

There is a small but useful improvement to the pre- 
vious sequence: instead of adding a non-selective hard 
pulse to convert the contra-axial xy-rotation into a selec- 
tive xy-rotation, we can apply a non-selective z-rotation, 
that acts equally on both spins, to reset one spin. 
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For the spin that we wish to return to its original state, 
the non-selective z-rotation cancels the z-rotation due to 
precession. Thus this spin does not, in effect, precess be- 
tween hard pulses and so is returned by the second hard 
pulse. Note that it is not actually necessary to perform 
the non-selective z-rotation; instead we can simply move 
into a new reference frame, effectively rotating the frame 
rather than the spins pjj. thus saving one pulse. Equiv- 
alently, we can pass the z rotation through the pulse se- 
quence, so that the selective Jump and Return sequence 
(equation [SJ becomes 
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where subscripts on pulses now indicate phase angles in 
the xy-plane. Each z rotation can be passed through 
successive pulses until it reaches the very end of the se- 
quence, at which point the final rotation can either be 
implemented using a composite z-pulse 'l8'| or, if the 
calculation ends with qubits in an eigenstate, simply be 
dropped. In some cases it is preferable to develop the 
pulse sequence as a series of independent modules, and 
in this case z rotations should be moved to the end of 
each module and then corrected using composite pulses. 

These techniques can be extended to provide a wide 
repertoire of single qubit gates. By the combination 
of contra-axial z-rotations (free precession) and non- 
selective z-rotations (frame rotations) we can perform 



selective z-rotations on either spin. In this way we can 
selectively rotate one spin about the z-axis to perform 
equally long xy-rotations on both spins, but around dif- 
ferent xy axes for each spin, 
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Using methods of this kind any desired combination of 
rotation and phase angles can be implemented. Alter- 
natively, such combinations can be constructed by using 
two selective pulses in succession. 



III. SOLVENT SUPPRESSION SEQUENCES 

The Jump and Return sequence is the simplest mem- 
ber of a large family of solvent suppression sequences [l^ 
developed for use in conventional NMR. Although Jump 
and Return does not excite signals exactly on resonance, 
this selectivity is not very robust, and signals only a small 
distance from resonance are significantly excited. A fam- 
ily of sequences, known as binomial solvent suppression 
sequences [T^ . have been developed which have much 
broader nulls in their excitation profiles. 

The best know of these sequence is usually called 1331 
and for a 90° rotation has the form 



(8) 



where the pulse lengths are chosen such that 8a = 90°, 
and r = 7r/2(5a; as before; the bars over pulse lengths in 
the name correspond to the changes in the pulse phases. 
Although this sequence is in some ways related to Jump 
and Return its basic operation differs in an important 
way: in Jump and Return sequences the target spin is 
fully excited by the first 90° pulse, while in binomial 
sequences the excitation builds up steadily throughout 
the pulse sequence. This leads to the use of a precession 
time 2t, compared with r in Jump and Return. 

As binomial sequences give more robust solvent sup- 
pression than Jump and Return, it might seem that they 
would also give more robust selective pulses, showing less 
dependence on small deviations from ideality in their im- 
plementation. There are, however, three reasons to doubt 
this surmise. Firstly, the binomial family is designed to 
give robust non-excitation of the solvent signal, not ro- 
bust excitation of the target spin. Secondly, the binomial 
family is designed solely for the purpose of excitation, and 
it is not clear how well it will behave as a general rotor 
|2Clj |. Finally the much longer evolution periods involved 
(6t for I33I, compared with r for a simple Jump and 
return sequence), mean that these sequences are likely to 
be vulnerable to errors in the value of (5a;, such as those 
arising from multiplet structure. 

We have investigated this question by modifying a 
range of binomial solvent suppression sequences to give 




FIG. 2: Propagator fidelity for a selective 90° pulse imple- 
mented using the Jump and Return approach of Eq. 2] as a 
function of the ratio g between the actual offset frequency 
and its nominal value. The solid line shows the fidelity with 
respect to a 90° pulse, while the dashed line shows the fi- 
delity against an identity operation. These curves are given 
by cos[7r(l =F5)/8] respectively. 



selective excitation sequences, and then calculating prop- 
agator fidelities as a function of small errors in buj. 
Although the detailed results are quite complex [23| , the 
overall result is simple: selective excitations based on 
Jump and Return sequences are far more robust than 
those based on more sophisticated solvent suppression 
sequences. For this reason we will only consider Jump 
and Return sequences for the remainder of this paper. A 
fidelity plot for one particular Jump and Return sequence 
is shown in Fig. [21 



IV. OFF-RESONANCE EFFECTS 

Although Jump and Return based sequences are gen- 
erally quite successful when implementing single qubit 
gates, there is a fundamental problem in their design. 
Since we use the difference in resonance frequencies to 
generate contra-axial rotations, it is obviously impossible 
for our hard pulses to be perfectly on-resonance with both 
transitions. In general we can parameterise off-resonance 
effects by the off- resonance fraction j^, / = Jw/wi, 
where (5w is the off-resonance frequency, and lo\ is the 
nutation rate induced by the RF pulse. The correspond- 
ing propagator for an x rotation is 



U = exp[-icji<(/2; + 

cos(C) - i! sin(C)/7 ~i sin(C)/7 

-isin(C)/7 cos(C) -ht/sin(C)/7 



(9) 



where <^ = ^ujit/2 and 7 = 
reduces to the simple form 



+ p. When / = this 



cos(a-'it/2) — i sin((jJit/2) 
-i sin(cjii/2) cos(wii/2) 



(10) 



as expected. When / ^ a pulse will induce a rotation 
around an axis tilted away from the xy-plane as shown in 
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FIG. 3: OfF-resonance effects on a 90° excitation pulse for 
off-resonance fractions / in the range to 1 in steps of 0.1. 



Fig. El However, the undesired effect of rotating around 
this tihed axis is a systematic error and can be compen- 
sated for by various means, including composite pulses. 

The use of composite pulses to correct for off-resonance 
errors has been widely studied in conventional NMR [23| . 
More recently, corrective composite pulses that act as 
general rotors, known as Class- A composite pulses p^ . 
have been developed for conventional NMR experiments 
psf and for NMR quantum computing . For a system 
with just two spins, with equal and opposite known reso- 
nance offsets, there is an analytical solution for a compos- 
ite pulse that perfectly corrects the off-resonance error. 
This ROTTEN (Resonance Offset Tailoring To Enhance 
Nutations) composite pulse j2Qj consists of a time sym- 
metric sequence of three hard pulses, that carries out the 
desired a;y-rotation on both spins. The fidelity of the de- 
sired rotation at the actual resonance offset of both spins 
is unity, since the rotten composite pulse sequence is 
tailored to perform perfectly for this resonance offset. 

Although it is possible to combine the three period 
ROTTEN sequence with the Jump and Return approach, 
a simpler sequence results if two of the RF pulses are 
converted to z-rotations. A helpful insight into how this 
new sequence works is provided by an existing result from 
conventional NMR. It is well known that a 90° pulse with 
a small off-resonance fraction / applied to an initial z- 
state (usually called an excitation pulse) is, to a good 
approximation, equivalent to an ideal 90° hard pulse fol- 
lowed by a period of free precession for a time equal to 
2/7r « 0.64 of the pulse length "2^. This can be seen in 
Fig. El where the lines corresponding to small values of / 
all terminate near the xy-plane, and are offset from the 
y axis by angles that are roughly proportional to the off- 
resonance fraction. The same conclusion can be reached 
by expanding Eq. Elto first order in /. 

This approach is not confined to excitation but can be 
generalized: a 90° pulse can be approximated as an ideal 
pulse surrounded by free precession periods of length 
equal to 2/7r times the pulse length. (The case of an 
excitation pulse, described above, is special: the initial 
precession period can be dropped as z-states are invariant 
under free precession, allowing a simpler approximation 



to be used.) Reversing this, it is clear that a 90° pulse 
applied to a spin with a small off-resonance fraction / 
would behave almost like an ideal pulse if it were pre- 
ceded and followed by periods of length equal to — 2/7r 
times the pulse length; methods for implementing nega- 
tive time periods are considered below. 

For larger values of / this approximation begins to 
break down in two ways. Firstly the phase angle is no 
longer linear in the offset frequency, and so it is not pos- 
sible to correct this using a fixed evolution time, and 
secondly the final state can lie some distance from the 
xy-plane. It is, however, possible to achieve an ideal 
rotation for any particular off-resonance fraction in the 
range |/| < 1 by placing negative z-rotations before and 
after a pulse with a nominal nutation angle greater than 
90°. The z-rotation sandwich 

0; C (11) 

(where 9' is the nutation angle of the central pulse in the 
absence of off-resonance effects) can be shown to perform 
an ideal 90a; rotation as long as the angles are chosen as 

<^' = -arcsin(/) 0' = ^^^^^^^±£1 . (12) 

Pulses with other phase angles can, of course, be achieved 
by changing the phase of the central pulse. For small off- 
resonance fractions Eg. 1121 can be expanded to first order 
in /, giving 0' « — / and 9' « 7r/2. Outside the range 
/I < 1 the angles 6' and (p' become imaginary, and this 
approach is no longer possible. 

This approach can also be used for two spins with equal 
and opposite resonance offsets; the nutation angle 9' will 
be the same for both spins, while the z-rotation angle 
will take equal and opposite values for the two spins; the 
two spins need a period of free precession given by 

t'^^X tgo (13) 

(where ^go is the length of the 90° pulse) to rotate in op- 
posite directions. For small values of /, t' = (— 2/7r)/:go 
as expected. The spin precession time is negative, but we 
can achieve this by subtracting the desired negative pre- 
cession time from any neighbouring positive precession 
period. If no such period exists then one can be created: 
the two spins will precess through angles of ±360° dur- 
ing a time 4r, and such time periods can be inserted with 
impunity. In effect this allows the spin to achieve a neg- 
ative precession angle by going the long way around the 
Bloch sphere. In fact it is not necessary to go all the way 
round, as a similar effect can be achieved by combining 
a precession period of length 2r, which gives rotations of 
±180°, with a simultaneous 180z rotation, so that one 
spin is left alone while the other rotates through 360°. 

The strategy of a z-rotation sandwich also works for a 
pulse with an arbitrary nutation angle, 9x. The required 
rotations are now given by 

(j)" = - arcsin(/ tan[6l/2]) (14) 
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and 

, _ arccos(-/^ + [1 + /^] cos 9) 



with first order approximations (/)" « — /tan(0/2) and 
9" « 9. Solutions exist for off-resonance fractions in the 
range |/| < cot(6'/2); note that as 9 increases, the range 
of off-resonance frequencies which can be tolerated de- 
creases. This is particularly obvious for a 180° rotation 
(a NOT gate), which can only be formed when / = 0, as 
any tilt in the rotation axis will make it impossible to 
rotate a state from the north pole of the Bloch sphere to 
the south pole. At the other extreme it is possible to use 
^-rotation sandwiches to perform accurate 5° rotations, 
even with off-resonance fractions greater than 20! 

A solution to this problem is to implement pulses with 
large rotation angles as two or more pulses with smaller 
rotation angles. For example it is possible to implement a 
180° pulse by combining two 90° pulses, which are them- 
selves implemented using z-rotation sandwiches. Thus 

r" 9'^ (4r + 2t") 9'^ r" (16) 

where r" = {2(f)" /it f) x tgg will implement a 29x pulse for 
two spins with equal and opposite off- resonance fractions. 
The relationship of r" to 9" and </>" is analogous to that 
between r', 9' and (/>' 

The z-rotations at either end of the sequence are back- 
ward precession times, but it may be possible to combine 
these with surrounding precession periods. The range of 
/ values which can be tolerated will be at least doubled 
by this process, and for large angles the improvement 
will be much greater. The process can, of course, be ex- 
tended indefinitely, dividing a pulse up into a very large 
number of smaller pulses; the final result is similar to 
a DANTE pulse 9], but for such large off-resonance fre- 
quencies it is probably more sensible to use conventional 
shaped pulses. 

For the special case of a 180° pulse (simultaneous not 
gates on both spins), the outer z-rotations can simply 
be dropped. In general a 180° rotation refocuses equal 
periods of precession on either side of the pulse, so we 
can add periods of precession to either side and it re- 
mains a 180° rotation. Thus we can "add" periods of the 
right length to exactly cancel the backward z-precession 
periods, and the sequence 

9', [At + 2r') 9', (17) 

will implement a 180^ pulse for two spins, and is valid 
for the range |/| < 1. 

This sequence can be shortened by replacing the cen- 
tral ±360° evolution by a ±180° evolution and a 180° 
frame rotation, to give the sequence 

9'^ 90^ (2r + 2t') 90^ 9'^ (18) 

where the 180° frame rotation has been split into two. 
For the permitted range of off-resonance frequencies 



(I/I < 1) the evolution period will always be greater 
than or equal to zero. The frame rotations can be passed 
through the pulses to the beginning and end of the se- 
quence, and can then be dropped. Thus the sequence 

9'_y (2r + 2r') 9'^ (19) 

will also implement a 180a; inversion pulse for two spins 
with I/I < 1. 

This approach can also be combined with conventional 
Jump and Return sequences, to perform more accurate 
selective or contra-axial rotations, by simply replacing 
each 90° pulse in a sequence with a z-rotation sandwich. 
As usual it is possible to combine the different periods of 
precession to give a simpler sequence. Thus using equa- 
tion and choosing 9 — 90° gives 

t' 9'_y (r/2 + 2t') 45, 9'y r' = [90, , 1] (20) 

where the central evolution periods have been combined, 
and the 45, rotation can be implemented using frame 
rotations. The precession times at either end of the se- 
quence are negative times, but as before it will be possible 
to combine these with surrounding precession periods or, 
failing this, to introduce additional precession periods of 
4t with which they can be combined. 

Once again it is instructive to consider the special 
case of a 180° contra-axial rotation, when the outer z- 
rotations can simply be dropped as before to get 

9'_yi2T + 2T')9'y (21) 

a sequence identical to that derived previously, equa- 
tion^] This makes sense: a 180° contra-axial rotation 
is indistinguishable from a 180° hard rotation, as when 
going half way round the Bloch sphere it does not matter 
which direction is taken. 



V. SUMMARY 

Jump and Return sequences provide a simple but ef- 
fective way of implementing frequency selective pulses 
in two spin systems. They have already proved useful in 
NMR quantum computation, and could also be applied in 
related systems where implementing conventional shaped 
pulses is impractical. The simplest approach is only suit- 
able when the difference in resonance frequencies is small 
compared with the nutation rate of the external field, 
but more complex approaches can be used to implement 
both selective and non-selective pulses in spin systems 
with larger off-resonance fractions. 

Given the power of existing shaped pulse approaches 
it is important to consider what advantages are offered 
by Jump and Return sequences, which can be summa- 
rized as simplicity and time-efficiency. Although modern 
NMR spectrometers are usually capable of implement- 
ing complex shaped pulses, older equipment has much 
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more limited capabilities, and Jump and Return pulses 
have been used to implement NMR quantum comput- 
ing on such devices \l2L . More importantly electron 
spin resonance (ESR) spectrometers are not yet capable 
of implementing complex shaped pulses, and so simpler 
approaches will be needed if current quantum computing 
experiments based on ESR ^2f\ and electron-nuclear dou- 
ble resonance (ENDOR) ,28|J are to be extended. Jump 
and Return sequences based on equation 0] are partic- 
ularly simple, requiring only pulses with simple phase 
angles and delays. Finally, Jump and Return sequences 



are also time-optimal in that they are the shortest pos- 
sible sequences capable of achieving selective excitation, 
thus minimising decoherence effects during pulses. 
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